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Abstract. The possibility of the crossover from the BCS pairing to the Bose-Einstein condensate 
(BEC) of diquarks with going down in density is discussed in the framework of in the Nambu 
Jona-Lasinio (NJL) model. We find that the quark matter at moderate density may be close to the 
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intermediate of the crossover, the precursory regime to the BEC phase. 
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^ INTRODUCTION 

It is now established that the ground state of QCD at extremely high baryon density is in 
the color- flavor locked superfluid [1] where quarks with all the three flavors participate 
in the Cooper pairing. The appearance of the superfluid CFL is completely due to the 
BCS mechanism saying that any arbitrary weak attraction between quarks leads to the 
tachyonic Cooperon on the Fermi surface at sufficiently low temperature. 

For the compact star phenomenology, however, one needs to consider other candidates 
with less symmetric pairing because the strange quark mass is not so much smaller 
compared to quark chemical potential which is at most of order 500 MeV even at the 
center of stars. The finite value of strange quark mass causes a stress on the CFL phase. 
This kinematic effect is known to bring about a rich variety of phases at moderate 
densities; see [2, 3] for the NJL model studies of QCD phase diagram. 

In addition to such kinematic effect, there is the other key ingredient which plays an 
important role at low densities; that is the strong coupling nature due to the asymptotic 
freedom of QCD. Going down in density, quark-gluon interaction becomes large, and 
this may lead a modification in superconducting phases. In fact, the ladder QCD calcu- 
lation of the coherence length indicates that the Cooper pair size decreases significantly 
toward low density and it could be of order of inter-quark spacing at /i = 500 MeV [4, 5] . 
This strongly suggests the crossover from the BCS-type weak coupling superconductiv- 
ity to the Bose-Einstein Condensate (BEC) of tightly bound quark pairs [6, 7, 8, 9, 10]. 
Another possible interesting phenomenon is the formation of pseudogap above the crit- 
ical temperature [11]. These two scenarios indicate the existence of non-trivial (non- 
Wigner) phase above the critical temperature at strong coupling. 

Major studies done in past mainly concentrate on the spectral analysis of collective 
modes. If the fluctuation is so large, however, there must be its feedback to the thermody- 
namics, the equation of state, for example. Such strong modification of thermodynamics 
due to fluctuations may bring about remarkable astrophysical consequences. 



In this short article, after briefly summarizing the application of the Nozieres Schmitt- 
Rink (NSR) theory [12] to the NJL model following our detailed analysis [6, 7], we dis- 
cuss the relativistic BCS-BEC crossover paying a particular attention to the significance 
of the fluctuation feedback to the quark matter thermodynamics. 

APPLICATION OF THE NOZIERES SCHMITT-RINK THEORY 
TO THE NAMBU JONA-LASINIO MODEL 

We here consider a general relativistic four-fermion model with a point attraction [6], 

Sf[\if,\jf\ = \\f{iB -m + ron) \j/ + G\j/UY 5 C\j/*\j/ T ir5C\j/, (1) 

where y, m denote the Dirac spinor field and its mass, /! is the chemical potential 
to adjust the asymmetry between particle and antiparticle, and G parameterizes the 
strength of attraction. The extension to the NJL model with isospin doublet, three 
colors is straightforward and the results will be given later. After introducing Hubbard- 
Stratonovich fields A(t, x) for iXj^ysCxj/, the fermion can be integrated out: 

Z = Z J @A3>A* exp (-Seff [A, A*] ) . (2) 

Zq = e^/^oO^) i s (-^g f ree f erm i on of the partition function, while S e ff[A,A*] 
is the effective action for the collective bosonic fields. According to [12], we in- 
clude the effect of fluctuation up to the second order in A. We have 5 e ff[A,A*] = 

TY.nl J^k)^ \-h ~~ Z(*'C0 M ,p)] I A(/ce»„, p) | 2 where 0) n is fermionic Matsubara frequency. 
Then the field A can be integrated out and the pressure of the system leads to 

p(^T)=p (^,T)+p Mc (^,T), (3) 

where p (n, T) = 2TY a =± I j^j In (l + e -\VF+"?-°n\/ T ) is the free fermion contri- 
bution, and p fluc corresponds to the fluctuation contribution defined by 

^ta(l-Gz M> r(iQ*,K)). (4) 

Here, Xp,T(i&N,K) = F.T.(T T [v/ r z'7 5 Cv/(T,x)] [y/ T iy 5 Cy ^(0,0)] f ) is the Cooperon at 
one loop, and £1^ is bosonic frequency. (See [6, 7, 1 1] for the explicit expression.) Note 
that, since we are approaching T c from above, the overall factor 1/2 is dropped because 
the phase and amplitude fluctuation contribute equally to the partition function. 

When the temperature (chemical potential) is decreased (increased) from the normal 
phase, the Cooperon becomes tachyonic at some critical point. The thermodynamic 
stability requires 1 — G^.r(0,K) > 0. Because the function takes minimum at K = 
for the system without density imbalance, the condition of criticality is simply 



l-G* Mc ,r c (0,0)=0, 



(5) 



which is nothing but the Thouless criterion. This condition generates a one-dimensional 
line in the (jU, T) -plane which we call the Thouless line. To see how large the fluctuation 
effect on the thermodynamics is, it is better to move on to the canonical ensemble. This 
corresponds to determine the value of }l c by means of 

-^(^T c ) + ^ r (^T c ) = — 2 , (6) 

where kf parameterizes the total density of the system. The first term is the free fermion 
contribution while the second term is the fluctuation contribution. We will see later that, 
the second term gives a significant contribution even in the BCS side. 

From above two basic equations, we numerically obtain T c and jl c as a function 
of density and coupling, i.e., T c (G,kp,A) and jj, c (G,kf,A) where A is an appropriate 
momentum cutoff. Before going into numerical computations, let us briefly summarize 
the effect of including N c = 3 colors and Nf = 2 flavors. If we assume the attraction in 
the isoscalar and color anti-triplet channel, the Thouless criterion Eq. (5) is not affected. 
On the other hand, the number condition Eq. (6) is modified as follows. 

N fN ^., T C ) + *fcl> Tc) = NfNc ± (7) 

k 3 

where kf is redefined so that each fermion species has a density when the interaction 
is turned off. We see that the fluctuation contribution is multiplied by a kinematic factor 

dg = Nc ^~ l \ representing the fact that the system has dg collective modes belonging 
to the antisymmetric representation of SU C (3). From the parametric dependence, we see 
that the fluctuation dominates the thermodynamics in the N c — > °° limit. 

In numerical calculations, we set |^ = 0.2. In addition, we mainly study the crossover 

detail for a relativity parameter, = 0.2 x , leaving to see its density /mass dependence 
later. Further we use the modified coupling g = ^ instead of bare coupling G 2 . 

In Fig. 1(a), we show T c , fi c as a function of g. When g <> —0.2, T c is well approxi- 
mated by the mean field result indicated by 7mf- In the mean field approximation (with- 
out PflucX Mmf — Ef by neglecting a tiny correction of order T 2 /Ep. As a consequence, 
Tuf is determined almost by the Thouless criterion. As g is increased and the unitary 
point g = is approached, the mean field result starts to deviate from real T c . This means 
that fluctuation contribution in the number equation Eq. (7) grows gradually and it can- 
not be ignored anymore. In fact, the quark density coming from fluctuation grows as 
seen in fig. (b). When the coupling exceeds g = 0.07, \x c gets lower than fermion mass m 



Note that even in the known strong coupling system, the nuclear matter consisting of neutron and proton, 
this parameter is of order ^ = i^(10~ 3 ) at the Mott transition point [13]. 

2 The regularized coupling Gr is introduced by — ^ = ^ — jo,o(2m,0), and Gr c stands for the critical 

regularized coupling for zero mass boson at vacuum, i.e., — ^ = jo,o(0,0) — jo,o(2?w,0). Gr is related 

to the scattering length a s by Gr = Then we see that the zero-binding bound state forms at g = 
(the unitary limit) and it becomes massless at g = 1 at vacuum. See [7] for the detail. 
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FIGURE 1. (a) T c and jj. c as a function of the modified coupling g. These are shown in the unit of the 
Fermi energy, Ep = \Jm 2 + k 2 F . For other lines, see text, (b) The quark number contents; the line denoted 
by Nq (A% uc ) corresponds to the ratio of the mean field (fluctuation) contribution to the total density. 



and the in-medium bound state appears accordingly. (See [6, 7] for the detailed spectral 
analysis.) Then the system goes into the BEC phase where T c is nearly saturated to a 
constant and 7mf completely fails. This saturating behaviour suggests that T c is deter- 
mined by the number equation Eq. (7) because it does not explicitly depend on G. In 
fact, the in-medium boson mass is shown to be twice of /i c (Mb — 2jU c ) at T = T c , and 
provided M B 3> T c (nonrelativistic), Eq. (7) can be approximated by 

NN 4 N C (N C -1) f dK 2 (3 .,M B TM 



The factor 2 in the integrand comes from the fact one diquark consists of two quarks. 

This gives an approximation T c ~ T^ c = ^r/i (^"C(|)) 1 m^- In fi g- ( a )' this formula 
with Mb = 2jU c is tried by the thin line, which agrees very well with real T c . The 
rightarrow indicates the nonrelativistic strong coupling limit of T c , evaluated by Tj^, 
with Mb = 2m. In the current framework, T c does not saturate to this value and slightly 
increases due to the binding effect. This is a residual relativistic effect that the binding 
energy of diquarks can become as large as the order of the constituent fermion mass. 

As g is increased and g = 1 is approached, starts to fail and the system eventually 
goes into the new regime, the relativistic BEC (RBEC) [14]. This is because the nonrel- 
ativistic approximation 2\x c ^> T c is no longer valid there due to the large binding effect. 
In fact, the diquark mass is smaller than temperature, 2ji c <, T c , for g J> 1, and therefore 
anti-diquarks contribute to the thermodynamics. By only taking the stable boson and 

rt 1 n? ~N~f 

antiboson contributions in Eq. (7), we get an approximation T c ~ T^ c = j^y j^jj—^ 
in the same way as [14]. This formula with Mb = 2\i c is tried by thin line in FIG. 1(a), 
which fairly agrees with the real T c . However, in contrast to the RBEC of the elementary 
boson [14], our composite boson system has fermionic degrees of freedom due to the 
competition between the internal energy and entropy [7]. For this reason, the agreement 
is not so good as that in the BEC region. FIG. 1(b) confirms this picture; quark density 




from the fermion sector again comes to play a major role in the deep RBEC region. 

Remarkably, in the (R)BEC region, there is a non-trivial phase above T c , the pre- 
formed boson phase, up to 7di ss . 7di ss , shown by the dashed line in FIG. 1(a), charac- 
terizes the ionization of diquarks. Needless to say, the system in the preformed boson 
phase differs much from a pure Fermi gas although the symmetry is restored. 

Interestingly enough, the definite crossover appears also in the dynamic equation for 
the pair excitaion [7]. By expanding S e ff up to the quartic order in A and performing the 
low energy/long wavelength expansion, we obtain the dynamic equation near T c 



d . 5F eff [A,A*] 
-id— A(f,x) — 



dt v ' ' 8A*(t,x) 



«o\^ + £v^ |A(f,x)| 2 



A(f,x), (9) 



where d (complex) and {ao,c,bo} (real) are the low energy coefficients. If we define 
the complex effective mass by M c g = its real part coming from the particle-hole 
asymmetry represents the propagating piece of the fluctuating pair field [16], while the 
imaginary part expresses its diffusive nature. The real and imaginary parts of M e g- as a 
function of g are depicted in FIG. 2(a). In the BCS regime, the pair mode is diffusive, 
but the magnitude of damping decreases significantly as the unitary point g = is 
approached. When the system goes into the BEC phase, the imaginary part vanishes 
and fluctuation becomes a pure propagating mode due to a bound state gap. However 
M e ff does not saturate to 2m in contrast to the nonrelativistic calculation [15]; it gets 
smaller towards the RBEC phase, which is also due to the relativistic binding effect. 

Fig. 2(b) shows how the crossover characteristics of T c is affected by the increase 
of relativity parameter This corresponds to decreasing m or, increasing density k F . 

From bottom to top, — = 0.2 x n with n = 1, 2, • ■ • , 12. Several notes are in order, 
(i) The BCS/BEC crossover point m — jl c indicated by the large point shifts to higher 
value of g. This is due to the Pauli-blocking effect which prevents the formation of in- 
medium bound state at high (small) density (mass). At the same time, the crossover 
characteristics of T c gets smeared, (ii) The universal thermodynamics of the unitarity 
point is absent in the relativistic system as noticed in [9]. T c explicitly depends on the 



additional parameter In fact, the pressure at zero temperature is deduced to take a 

form p = j!^^f(-£^f, 7^"> ^) w i m / denoting a dimensionless function; we observe 
only in the m — > °o limit the nonrelativistic universal behaviour at = is recovered. 

Let us finally discuss in which regime the actual quark matter does exist. The modified 
coupling g is a function of A and m. It is easy to see g — > 1 with m — > °°. Also g — > — °° 
when m — > with a natural assumption that a bare coupling GA 2 is less than n 2 , the 
critical value for the dynamical Majorana mass generation at vacuum. As there should 
be the point g = in between, we conclude that for any fixed value for G < % 2 , the 
system is BCS-like for m — > and it is BEC-like at sufficiently large m. If we fix G to 
the usually adopted value, 3/4 of an appropriate scaler coupling [2, 3], g = corresponds 
to m/A = 0.53. This is much larger mass compared to the current quark masses in 
agreement with the recently appeared paper [16]. We conclude that somewhat exotic 
conditions must be satisfied to have the diquark BEC in QCD phase diagram; (i) the 
diquark coupling is stronger than expected, and/or (ii) in-medium quark mass is larger 
than its perturbative estimate. Interestingly, two lattice studies, one about qq interaction 
[17] and the other for in-medium quasiquark mass [18] are encouraging. Also it is worth 
noting that even if the BEC phase cannot be reached, the fluctuation feedback to quark 
matter thermodynamics will be significant. In fact, it can be neglected only in the weak 
coupling limit as noted in [10]. If so it may bring about a remarkable modification of the 
structure of possible quark, or hybrid stars. Exploring possible BCS/BEC crossover in 
quark matter with more realistic situations taken into account, as well as looking for its 
astrophysical consequences clearly deserves further investigations. 

This work was supported by a Grant- in- Aid for the 21st Century COE "Center for 
Diversity and Universality in Physics". The numerical calculations were carried out on 
Altix3700 at YITP in Kyoto University. 

REFERENCES 

1. K. Rajagopal and F. Wilczek, arXiv:hep-ph/001 1333; M. G. Alford, Ann. Rev. Nucl. Part. Sci. 51, 
131 (2001); T. Schiifer, arXiv:hep-ph/030428 1 ; D. H. Rischke, Prog. Part. Nucl. Phys. 52, 197 (2004). 

2. H. Abuki, M. Kitazawa and T. Kunihiro, Phys. Lett. B 615, 102 (2005) [arXiv:hep-ph/0412382]; 
H. Abuki and T. Kunihiro, Nucl. Phys. A 768, 118 (2006) [arXiv:hep-ph/0509172]. 

3. S. B. Ruester et al, Phys. Rev. D 72, 034004 (2005); D. Blaschke et al, Phys. Rev. D 72, 065020 
(2005) 

4. M. Matsuzaki, Phys. Rev. D 62, 017501 (2000) [arXiv:hep-ph/9910541]. 

5. H. Abuki, T. Hatsuda and K. Itakura, Phys. Rev. D 65, 074014 (2002) [arXiv:hep-ph/0109013]; 
K. Itakura, Nucl. Phys. A 715, 859 (2003). 

6. Y. Nishida and H. Abuki, Phys. Rev. D 72, 096004 (2005) [arXiv:hep-ph/0504083] . 

7. H. Abuki, Nucl. Phys. A 791, 1 17 (2007) [arXiv:hep-ph/0605081]. 

8. K. Nawa, E. Nakano and H. Yabu, Phys. Rev. D 74, 034017 (2006); J. Deng, A. Schmitt and Q. Wang, 
arXiv:nucl-th/061 1097; A. H. Rezaeian and H. J. Pirner, Nucl. Phys. A 779, 197 (2006). 

9. L. He and P. Zhuang, Phys. Rev. D 75, 096003 (2007) [arXiv:hep-ph/0703042]. 

10. L. He and P. Zhuang, arXiv:0705. 1634 [hep-ph]. 

11. M. Kitazawa et al, Phys. Rev. D 65, 091504 (2002); Phys. Rev. D 70, 056003 (2004); Prog. Theor. 
Phys. 114, 117(2005). 

12. P. Nozieres and S. Schmitt-Rink, J. Low. Temp. Phys. 59 195 (1985). 

13. U. Lombardo etal, Phys. Rev. C 64, 064314 (2001). 



14. H.E. Haber and H.A. Weldon, Phys. Rev. Lett. 46, 1497 (1981); J.I. Kapusta, Phys. Rev. D 24, 426 
(1981). 

15. R. Haussmann, Phys. Rev. B 49, 12975 (1994). 

16. M. Kitazawa, D. H. Rischke and I. A. Shovkovy, arXiv:0707.3966 [nucl-th]; arXiv:0709.2235 [hep- 
ph]. 

17. A. Nakamura and T. Saito, Prog. Theor. Phys. 112, 183 (2004) [arXiv:hep-lat/0406038]. 

18. P. Petreczky et ai, Nucl. Phys. Proc. Suppl. 106, 513 (2002). 



